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Abstract 

A cut on the maximum lifetime in a lifetime fit does not only reduce 
the number of events, but it also, in some circumstances dramatically, 
decreases the statistical significance of each event. The upper impact 
parameter cut in the hadronic B trigger at CDF [l] [2] [3], which is due 
to technical limitations, has the same effect. In this note we describe and 
quantify the consequences of such a cut on lifetime measurements. We 
find that even moderate upper lifetime cuts, leaving event numbers nearly 
unchanged, can dramatically increase the statistical uncertainty of the fit 
result. 
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Figure 1: Given the 3- momenta of all particles in the decay, the cut on the 
Impact parameter of the decay products translates directly into a cut on the 
dccaylength and hence on the lifetime of the primary particle. For clarity, the 
figure only illustrates the effect of an impact parameter cut on one of the decay 
products (the one going straight upwards). 
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1 Introduction 

In this note we discuss the impact of an upper lifetime cut on the precision of a 
lifetime measurement. We find that even an upper lifetime cut that loses only 
a small fraction of the data can have dramatic consequences on the precision of 
the lifetime fit, due to a loss of statistical power of those events that pass the cut. 
The small loss of events due to a moderate upper lifetime cut is accompanied 
by a large loss of information, because not only a few events outside the allowed 
time window are lost, but also the information that there were only a few. 
This can have dramatic effects on the precision of the measurement. As shown 
below, an upper lifetime cut that loses 14 % of the data reduces the statistical 
significance per event by 72 %, so the combined effect on the statistical precision 
of the lifetime measurement is equivalent to losing 76 % of the data. 

Such a cut on the maximum lifetime is for example implicitly applied in the 
hadronic trigger sample at CDF pQ, where the trigger requires two tracks with 
a minimum impact parameter of 100 — 120 /im (depending on the exact config- 
uration) and a maximum impact parameter of 1000 /^m [5] [3]. These impact 
parameter cuts translate into upper and lower lifetime cuts, which differ event 
by event. This is illustrated in Figure [1] 



2 Likelihood with lifetime cuts 

We can write the probability to find an event with decay time t, given that it 
passed the trigger cuts, as: 
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where we ignore the effect of measurement errors. The total log-likelihood func- 
tion for a set of N "ideal" two-body decays (no measurement uncertainties, 
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background, etc) is given by: 



log£ = -N log (r) 
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where the index i labels the event, each of which has its measured decay time 
ti and minimum and maximum lifetime cuts of t m i n i and t max i- 

Note that the only difference compared to the the likelihood function without 
lifetime or impact parameter cuts is the term: 
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Including Gaussian event-by-event measurement errors, the PDF for an event 
with lifetime cuts t m i n and t max , and the measured lifetime to, is given by: 
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3 Calculating the Uncertainty of the Fit Result 

The variance on the fit result can be estimated as the negative inverse of the 2 nd 
derivative of the likelihood function, evaluated at the lifetime t that maximises 
the likelihood: 

s-v-J*m (6 ) 
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This formula assumes that near its minimum, the Likelihood curve can be ap- 
proximated by a Gaussian. While this approximation is not perfect for lifetime 
fits, that have asymmetric errors, it is still very good as long as the event samples 
are sufficiently large, as shown in Section [5] 

For simplicity, we ignore event-by-event measurement errors in our subsequent 
calculations. We will show in Section [5] that, for the purpose of estimating the 
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error on the fit result, this provides a very good approximation for the case of B 
lifetimes measured at CDF where the event-by-event measurement uncertainties 
are much smaller than the lifetime to be measured. The 1 st derivative of the 
likelihood function defined in Equation [21 is: 
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where we introduced Aij = i max $ — t m ; n , , the width of the time interval to which 
the i — th event is confined due to impact parameter, decay distance, or direct 
lifetime cuts. The 2 nd derivative is: 
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where the angle bracket indicate taking the mean of the expression inside over 
all events. At the value of r that maximises the likelihood, the first term of 
Equation [HI vanishes, and the variance is given by: 
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Note that the lower lifetime cut by itself does not have any impact on the 
statistical precision (apart from changing the number of events), it is the width 
of the time interval defined by the cuts, that matters; it is therefore the presence 
of an upper lifetime or impact parameter cut that affects the statistical precision 
per event. 



4 Statistical Power Per Event 

The right hand side of Equation [5] can be separated into two factors: 

• The variance in the absence of any upper lifetime cut: 
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• The correction factor due to the upper lifetime cut: 
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Figure 2: Statistical Power per Event as a function of At/r. 




V is defined as 



V = l-( 
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where At is the width of the 
time window defined by the 
lifetime cuts, r is the lifetime 
to be measured. 
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So the change in statistical precision per event due to an upper lifetime cut is 
accounted for by making the following replacement for TV: 

It makes therefore sense to define the statistical power per event, V, as 

^-(ys^ 1 >• <i3) 

V is 1 for events without upper lifetime cuts, and < 1 otherwise. It is defined 
such that TV events with an upper lifetime cut (where TV is the number of 
events after the cut has been applied) arc statistically equivalent to TV ■ V events 
without an upper lifetime cut. Figure [5] shows the statistical power per event as 
a function of At/r, the time interval defined by the cuts divided by the lifetime. 
Upper lifetime cuts that seem harmless at first sight, because they have a small 
effect on the number of events, can have a dramatic impact on the statistical 
precision of the lifetime measurement due the reduction in statistical power per 
event. For example an upper lifetime cut leaving a time interval that is twice as 
wide as the mean lifetime to be measured {At/r = 2) retains (1 — e~ 2 ) = 86% of 
the events, but the statistical power per event is reduced to 28%. The combined 
effect is equivalent to losing 100% — 86% • 28% = 76% of the unbiased sample 
before the cut, rather than the naively expected 14%. 

The CDF hadronic B trigger requires two tracks with impact parameters be- 
tween 100 [im and 1000 (im, which translate to different upper and lower lifetime 
cuts for each event, typically yielding a At between 1 and 3 times the B lifetime 
(this is an approximate number from studies of B° — ► irir decays) . So each event 
in that sample is, for the purpose of lifetime measurements, only worth about 
30% of an unbiased event. Note that this is true for lifetime measurements, only, 
and not for asymmetries or oscillation measurements, where it is the oscillation 
period that determines the scale At needs to be compared to, rather than the 
mean lifetime. 
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Figure 3: Fit to ~ 13k MC-generated signal events. The shaping of the distri- 
bution due to the trigger manifests itself as a clear deviation from a straight line 
in this logarithmic plot. The fit describes the data well with a x 2 /dof = 0.8. 
The result of cr = 459.11^ fj,m is in good agreement with the input value of 
462 /mi. The numerical error estimate of ^y'f/Am, calculated using the MINOS 
algorithm in MINUIT, agrees well with analytical value of ±7.2 /mi from Equa- 
tion [5] d erived here. 




lifetime xc 



5 Monte Carlo Studies 

In order test how good the error estimate in Equation [5] is, given that this 
simplified formula ignores measurement errors and assumes the likelihood to 
be Gaussian near its minimum, it is compared to the error estimate from a 
MINUIT fit using a likelihood function that includes event by event errors. The 
errors are calculated using the MINOS algorithm within MINUIT. The fit is 
performed on simulated B° — > tttt events at CDF. The Monte Carlo simulation 
includes a detailed description of the CDF detector, including the hadronic B 
trigger, that requires two tracks with impact parameters between 100 /mi and 
1000 /im. This requirement translates into upper and lower lifetime cuts, which 
differ event by event. The likelihood function used to fit the simulated data 
includes the trigger effects, and the event-by-event uncertainties in the lifetime 
measurement, as given in Equation |U It is in the calculation of the acceptance 
limits where further detector effects arc taken into account, in particular the 
difference between the fast-reconstructed online impact parameters used by the 
trigger, and the more precise offline reconstruction used in the actual lifetime 
measurement. The result of the fit to 12,711 signal events with an input 
lifetime of cr = 462 /mi, is ct = 459.1 /im. The agreement between the fit 
and the MC data is excellent, as can be seen in Figure [3] 

The average statistical power per event, calculated using Equation I13[ is 31 %. 
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Using this in Equation yields an error estimate, ignoring measurement errors, 
of 7.2 /im, in good agreement with the numerical error estimate from the full 
likelihood. The solid line in Figure |4] represents the — 2 log(likclihood) curve for 
this fit, normalised to — 21og£ = at its minimum. The broken line represents 
the parabolic (— 2 log (Gaussian)) approximation implicitly made in Equation^ 
The agreement between the parabolic approximation and the actual likelihood 
for our data sample of 13k events is very good within the range relevant for 
calculating la errors, i.e., for this normalisation, up to (— 21og£) = 1. The 
errors estimated using the full likelihood are +7.3/im, — 7.1/xm, in excellent 
agreement with the parabolic approximation of ±7.2/ito. The difference between 
the two curves increases for larger (— 21og£) values, as the asymmetry in the 
lifetime errors becomes more pronounced, but it remains reasonable over the 
entire range of the plot, up to (— 21og£) r* 40 and beyond. This implies that 
the parabolic approximation remains reasonable even for much smaller data 
samples containing only a few hundred events. 

This gives us confidence that the formula derived in Equation is correct, 
and that, for the purpose of estimating the statistical uncertainty of a lifetime 
measurement, the approximations we made are justified for B hadron lifetime 
measurements at CDF, where event samples are usually large and typical event- 
by-event lifetime errors are about 60 fs, small compared to B hadron lifetimes 
of about ~ 1.5 ps. 
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6 Summary 



We quantified the statistical effect of upper lifetime cuts in lifetime measure- 
ments for the simplified case that the event-by-event lifetime errors are small 
compared to the lifetime to be measured. 

We found that the effect of an upper lifetime cut is generally much more dra- 
matic than the mere loss of events would suggest. The greatest impact of such 
a cut is a reduction in the statistical significance of each event for the purpose 
of lifetime measurements. For example an upper lifetime cut at twice the mean 
lifetime to be measured loses only 14% of events, but the statistical power per 
event is reduced by 72%. The combined effect is equivalent to a reduction in 
sample size by a factor 4, thus doubling the statistical error. 

We verified our calculation using simulated B° — > tttt events at CDF. In the 
trigger scenario used as an example, the statistical power per event is reduced 
by 69% due to the upper impact parameter cuts applied by the trigger. 
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